A Statistical Mechanics Approach to the Inherent States of Granular Media 
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We consider a Statistical Mechanics approach to granular systems by following the original ideas 
developed by Edwards. We use the concept of "inherent states", defined as the stable configurations 
in the potential energy landscape, introduced in the context of glasses. Under simplifying assump- 
tions, the equilibrium inherent states can be characterized by a configurational temperature, 1//3. 
We link f3 to Edwards' compactivity and address the problem of its experimental measure. We also 
discuss the possibility to describe the time dependent distribution probability in the inherent states 

■ with an appropriate master equation. 
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The possibility to describe general features in the physics of granular media with the concepts of Statistical Me- 
chanics was suggested few years ago by Edwards 0. As much as systems of standard Statistical Mechanics, each 
macroscopic state of a granular medium corresponds to a huge number of microstates. Furthermore, these systems 
show very general reproducible macroscopic behaviors characterized by a few control parameters In granular 
, media, however, in absence of some external driving the microstates are "frozen", i.e., they don't evolve in time as 

■ thermal energy is negligible compared to gravitational energy. Therefore the external thermal bath can be considered 
at zero temperature. Thermal motion of grains can be replaced, instead, by agitation induced by shaking or other 

■ form of driving, which may generate a dynamics among "frozen" microstates. In facts, granular systems are just 
one example of a broad category of materials which can be found in "frozen states" . An other interesting example 
is the state of a supercooled liquid which has been cooled at zero temperature and falls in one of the microstates 
corresponding to the minima of the potential energy, called "inherent structures" || . Actually, the analogy between 

(-~) slow dynamics in granular media and glassy behaviors of supercooled liquids has been proposed in various forms 

Eli- 

In this Letter we will further extend this analogy to develop a statistical mechanics approach to granular materials, 
following Edwards pioneering ideas Jlpj ] . For simplicity we consider here a granular material as made of particles which 
may have a distribution of shapes and interacting via hard core potentials, so that they are not allowed to overlap. 
Due to gravitational energy, each particle configuration {r;} corresponds to a potential energy U{ri}. The stable 
packing configurations are the minima or saddle points of the potential energy or, more generally, all the states which, 
due for instance to friction, are mechanically stable. These states are the frozen configurations which by analogy 
with the glass terminology we call inherent states. If there is no ambiguity throughout the paper we call them simply 
Q states. We label each state by an index s and call E s the corresponding energy. 

The dynamics from one state to another can be induced by some form of driving such as, for instance, by sequences 
of shakes or "taps" of given amplitude as in typical experimental realizations (see references [||^]). The kinetic 
• i-H . energy driven in the packs by the "shakes" is totally dissipated when the shake is over and the system is almost 
' instantaneously "frozen" in one of the inherent states. The actual dynamics depends on several parameters such as 
the tap amplitude r = a/g (a is the shake peak acceleration and g the gravity acceleration), frequency and others 
jtjj . We say that the system is at stationarity if its macroscopic properties don't change any longer with respect to 
the given dynamics allowing transitions from one microstate to another. We refer to these states as stationary states, 
or quasi-stationary if the macroscopic properties change very slowly. 

A very important issue is to individuate the states distribution, namely what is the probability to find the system 
in a given inherent state s. We can expect that under stationary condition Jl^| the system has been "randomized" 
enough, and therefore, following essentially Edwards original ideas, we make the assumption that such a distribution 
is given by a maximization of the entropy under the condition that the average energy is fixed. 

More precisely we consider a statistical ensemble of equivalent systems all prepared in the same way. We indicate 
with {E r } the energies of the accessible inherent microscopic states of each system and with n r the number of systems 
with energy equal to E r . The average energy per system is thus E = ^ r P r E r where P r — n r /N is the probability 
to have a system in the inherent state r. 

The configurational entropy can be defined as S — — J2 r Pr m Pr ■ We assume that the stationary distribution is 
given by the maximal entropy under the constraint which fixes the average energy E. This requirement leads to the 
Gibbs distribution function: 

e -0E r 

Pr = — (1) 
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where the partition function Z = ^ r e~@ Er is a normalization factor and j3 a Lagrange multiplier determined by the 
constraint on the energy. 

As in standard Statistical Mechanics, we can show that in the thermodynamic limit S is the logarithm of the number 
of microscopic inherent states, f2, corresponding to the macroscopic energy E. Actually, we have that In Z = S — (3E. 
On the other hand Z can be written as Z = J2 r e~^ E - = ^(^)e" /3£ and using the steepest descent method, in 
the thermodynamic limit one obtains: InZ = InQ(E) — (3E. Thus, we have that: S = In 0(1?) and 

» = §§• 

We call = T con f the configurational temperature. We note that beside the configurational temperature T con f there 
is a second "temperature" related to the grains kinetic energy, the external effective bath temperature, Tbath, which 
is zero. While the fast degree of freedom, related to the momentum distribution, equilibrates with the temperature 
of the external bath resulting into a zero kinetic energy, the non zero configurational temperature characterizes 
the equilibrium distribution among the inherent states. Note that two granular systems in "thermal contact" have 
always the same Tbath = while the configurational temperature T con f may be different. However the configurational 
temperature characterizes the equilibrium states of the system and depends only on the average energy of the inherent 
states and not on the particular dynamics used. T con f enters as a parameter in the equilibrium distribution allowing to 
substitute, as in usual Statistical Mechanics, time averages, which are theoretically difficult to evaluate, with ensemble 
averages. 

The existence of two effective temperatures is quite familiar in the context of glassy systems []l3| , p~4| . If one quenches 
these systems in their glassy phase at a low bath temperature, Tbath, after a relative short transient the fast degree 
of freedom come quickly in equilibrium with Tbath, while the slow degree of freedom related to configurational rear- 
rangement fall out of equilibrium and are characterized by an internal temperature Ti nt |l4|]l5| ]. Interesting enough 
when Tbath — 0, Ti nt coincides with the configurational temperature introduced in eq.(Q): Ti nt — T con f CRjlq| . 

There are however some differences between glassy systems quenched at Tbath almost equal to zero and granular 
materials. In the former the energy depends on the equilibrium temperature before quenching and the dynamics 
among the inherent states is due to an activated thermal process at extremely low temperature. In granular media 
the energy depends on the shaking amplitude T and the dynamics in general proceeds with a tapping process with 
the same fixed amplitude. It would be interesting to study glassy materials using the dynamics of granular materials 
and vice- versa. For example one could introduce in glass formers a "tap dynamics" , in which the temperature Tbath 
is cyclically varied between zero and a finite value Tr |l7j . On the contrary in granular material one might first reach 
a given value of the energy by shaking the system with an amplitude T and then evolve the system by gently shaking 
with very low value of T. 

We now show that in the particular case in which the particles density p is constant, [3 can be easily related to the 
"compactivity" , introduced by Edwards in his seminal papers [0 . An "inherent state" can be characterized by the 
particle density function p{r). When the system has a constant density profile, i.e., p(r) = po in the packs, a one to 
one correspondence exists between energy and volume. Let's consider for sake of simplicity only hard core and gravity 
interactions. Actually, if A is the horizontal area of the container and h the packs height (its volume is V = hA), the 
system energy is: E = rngp^Ah 2 /2 = VV/2, where m is the grain mass and the characteristic pressure V = mgpoh is 
the weight of the sample per unit surface. Then we have: 

dE "P dV "p [ ) 

The quantity X^ 1 — ^ is Edwards' "compactivity" 0, which in this case is simply proportional to (3. 

To see how the formalism can work we consider a simple example where p is not a constant. As stated, the inherent 
states of a granular system can be characterized by their coarse-grained density spatial distribution, p(r), and the 
partition function of the system may be written as: Z = J Vp e _/3 ^ p ], where J-[p] — E[p] — f3~ 1 S[p] is the free energy 
functional for the system in its inherent states and the integral is over all inherent states. In the simplest cases the 
gravitational energy E[p] is given by: E g [p] = mg J v dr zp{f) (z is the component of r along gravity acceleration, g). 
S can be approximated by the entropy of a lattice gas of non overlapping grains with : 
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Here po is the maximal accessible density. The minimum of the free energy gives the following Fcrmi-Dirac equilibrium 
density profile: 

p(f) = PFD (z) = emgp Jl x) - 1 (5) 
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where A is a chemical potential which fixes the total number of particles. This result is in agreement with experimental 
HID and other theoretical findings. 

The configurational temperature can be evaluated by introducing a small perturbation jU where 7 is a small 
parameter and U is an observable. If U r is the value of U in the inherent state r the energy of the inherent state 
changes from E r to E r + jU r . Using the fluctuation dissipation relation 



d(U) 



(h. 



= P((U 2 ) - (U) 2 ) (6) 
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The ratio of the response function to the fluctuation of U gives (3. It is possible that such response function cannot 
be measured experimentally. Another way to link fj to observable quantities easier to measure, is the following 
From eq.(Q) using standard statistical mechanics immediately follows 

= (E 2 ) (E) 2 (7) 

which relates the fluctuation of the energy to the energy derivative. The average is over the equilibrium distribution 
of eq.(m). Both the energy and its fluctuation, C = (E 2 ) — (E) 2 , can be experimentally measured as function of T. 
Using -gpL — from eq. ([?]) we obtain 

where Tq is a reference point which could be adequately chosen. Therefore a part from an additive constant in 
principle it is possible to obtain /3 as function of T. This procedure was experimentally proposed in ref. We 
note that this approaches predicts an universal value of /?, independent on the particular control parameter T used in 
different experiments p"T[ . 

So far for simplicity we have considered the case in which the distribution is stationary. From the experiments Q] 
and from granular models p] we know that the dynamics often enters a quasi stationary regime in which average 
macroscopic quantities, like the energy or density, changes extremely slowly in time. In these situations we expect the 
distribution be quasi stationary as in glass models Jl5|, namely still given by (0) where (3 is slowly dependent on time. 
Here, (3 can also be defined from generalized dynamical fluctuation dissipation relations (FDR) extended to granular 
media j^o|. In analogy to eq.(^), in the linear response theory the FDR relates the system response to a perturbation 
to its time dependent correlation function in the unperturbed state |2l| , |22fl . 

The distribution given by eq.(|l|) is the equilibrium distribution. When the system is not in equilibrium the prob- 
ability P r will depend on the discrete time t n = nr where n is the n-th tap and r the duration of a single tap. In 
principle it is possible to write a master equation for such time dependent distribution: 

dP 

7 n L = -J2PrW rs + Y,P s W sr (9) 
OVfi 

s r 

where W rs are the transition probabilities from the inherent state r and the inherent state s. If the system has 
to approach the equilibrium distribution in the inherent states space we can assume that detailed balance must be 
satisfied: 



W r8 /W sr =exp[-/3(E S -E r )] . (10) 

Now we want to discuss how the statistical mechanics formalism discussed in this paper can work. In a realistic 
model for granular medium, in order to reproduce a tap process and then calculate the macroscopic quantities, one 
should consider the detailed dynamical, including the energy dissipated during the collision and then calculate the 
time average. The formalism developed here suggest to calculate the macroscopic quantities when the system has 
reached a stationary distribution, by taking ensemble averages instead of time averages. Exchanging time averages 
with ensemble averages is in fact one of the powerful feature of statistical mechanics. To construct the statistical 
ensemble one does not need to follow necessarily the exact dynamics but any dynamics which leads to stationarity. 
In principle one could evaluate P r (t n ) directly from (|l0|). Another possibility is to construct a dynamics which allows 
to visit the inherent states and mimics the tapping process. For instance, at the simplest level, a single tap can 
be schematically reproduced by introducing an effective time dependent bath temperature, Tt> a th(t), which is a step 
function which assumes two values: Tf, a th{t) = Tr > when the vibration is on (and the grains average kinetic 
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energy is larger than zero) and Tb at h{t) = when the vibration is over (and the kinetic energy is zero). A "tapping 
sequence" is the cyclical repetition of the above schematic tap [@J|]. This mimics the shaking process where Tr 
plays the role of the shaking amplitude T in the experiments, and allows to explore the inherent states visited when 
Ttath = 0. In general the parameter Tr is different from the temperature T con f previously introduced to describe 
stationary inherent states: Tr is related to the agitation of grains induced by the external drive, while T con f is the 
configurational temperature which characterizes the "frozen" inherent states. These two quantities are linked through 
the functional relation (||) where T is replaced by Tr- 

During the shaking process, when the kinetic energy is different from zero, the state k of the system is not necessarily 
a mechanical stable configuration. Say Vk(t) the probability that the system is in one of such states at time t. After 
each tap the system goes into one of the inherent sates, say r. Thus the probability, P r (t n ), to be in the inherent 
state r after the n-th tap is related to the the probability V r (t) by the following relation: 

Pr(t n ) = VrijlT) (11) 

So, within the above models it is possible to evaluate the time dependent distribution function for such particular 
dynamics and the corresponding energy. At stationarity say E the value of the average energy, the asymptotic 
distribution is given by eq.(|^) with (3 corresponding to the energy E. The above simplified approach to the dynamics 
of granular media has been recently also discussed in Ref. Q , with an interesting detailed discussion. 

Like in statistical mechanics the basic postulate that time averages can be substituted by Gibbs ensemble averages 
can be considered as a working hypothesis which cannot be proven but can be accepted only a posteriori if the results 
which follow are correct. 

There are many lattice models for granular media which are based on the above tap dynamics. They are 

based on an Hamiltonian formalism which takes into account gravity and hard core repulsions which may generate 
geometric frustration by not allowing two grains to overlap (for a review see ||). Although the models are rather 
crude, they have given general results in excellent agreement with experimental findings ranging from logarithmic 
compaction to segregation, density fluctuations, "memory" and "aging" effects, Fermi-Dirac density profiles, non 
Gaussian density distributions and others (for a review see 0). As an example, we show in Fig.^Jthe average density 
of the inherent states, p, in one of these models, the Tetris |J, as a function of Tr. The system is "shaken" with a 
given Tr for a long time and, after switching off the shaking, the density of the obtained "frozen" state is recorded. 

It is very interesting that the results of Fig.|| reproduce qualitatively the reversible line in granular compaction 
|Q| and closely remember the energy of inherent structures as a function of their equilibration temperature before 
quenching ||. 

In conclusion, in this paper we have elaborated the original Edwards' statistical mechanics approach to granular 
media to describe the quasi-static distribution of the system in their state at rest that we have called inherent states, 
in analogy with glassy systems. We have distinguished two different temperatures the zero external bath temperature 
and the internal configuration temperature related to Edwards' compactivity, which can be experimentally measured 
and linked to the driving T. We have postulated a Master equation, eq.(^|), for the evolution of granular packs among 
its inherent states. A straightforward application of these ideas reproduces the Fermi-Dirac density distribution found 
exp er iment ally. 

Models from standard statistical mechanics, which can be well understood in the present unifying formalism, are 
able to describe a variety of properties of granular media (l],|| @,||,[l9 20 2^j2^] . 
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FIG. 1. The average density of the inherent states, p, of a model of granular media (Tetris) as a function of the adimensional 
shaking amplitude Tr. The asymptote (the horizontal dashed line) corresponds to the minimal stability density p m . 
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